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Abstract 

A damage-plasticity constitutive model for the description of fracture in plain concrete is presented. Two 
approaches, the local model comprising the adjustment of the softening modulus and the nonlocal model 
based on spatial averaging of history variables, are applied to the analysis of a concrete bar subjected to 
uniaxial tension and to a three-point bending test. The influence of mesh size and the decomposition into 
damage and plasticity components are discussed. It is shown that for the two examples studied, both 
approaches result in mesh-independent results. However, the nonlocal model, which relies on spatial 
averaging of history variables, exhibits sensitivity with respect to boundary conditions, which requires 
further studies. 
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1 Introduction 

Concrete is a strongly heterogeneous and random material, which in the hardened state consist mainly 
of aggregates, cement paste, air voids and water. The heterogeneities range from several nano-meters 
in the structure of the cement paste to a couple of centimeters for aggregates. Macroscopic fracture in 
this heterogeneous material is characterised by a quasi-brittle response. Displacement controlled tests of 
the mechanical response of concrete in uniaxial tension, for instance, exhibit softening, which is defined 
in this context as a gradual decrease of the load with increasing displacements. This softening response 
is caused by the development of a fracture process zone, involving mechanisms such as branching and 
bridging of micro-cracks. The fracture process zone evolves from distributed and disconnected micro- 
cracks to a strongly localised tortuous macroscopic crack, whereby the tortuosity is governed by the 
presence of aggregates, which represent the dominant heterogeneity of the material. 

For the analysis of concrete structures, it is important to include the quasi-brittle softening response of 
concrete, since it strongly influences the load capacity of structures and is the source of a particular size 
effect on nominal strength, which neither follows the strength theory nor linear clastic fracture mechanics, 
see [l|. However, the heterogeneities in concrete described above, are usually significantly smaller than 
the size of structures analysed, which makes it often not practical to model the heterogeneities explicitly. 
Instead, concrete is usually assumed to be homogeneous on the scale studied and the formation of cracks 
is assumed to be an averaged process. For many heterogeneous materials, the averaged response can be 
obtained from spatial averaging of the material response on lower scales. However, for localised fracture, 
as in the case of concrete, spatial averaging alone is statistically insufficient, since energy is dissipated 
in localised zones. Thus, additional to spatial averaging, averaging of many individual fracture process 
zones has to be performed. This average of tortuous crack patterns results in a spatially distributed 
fracture process zone, which is governed by the size of the maximum heterogeneity. The fracture process 
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zone is denned as the average of the zone in which energy is dissipated during the fracture process. At a 
late stage of the fracture process, the fracture process zone in a single random analysis is fully localised 
in a single crack. The deterministic nonlocal model represents the average of many random analyses. 
Therefore, the average of fracture process zones at this stage of cracking is not localised but distributed, 
with its size determined by the tortuosity of crack paths. 

Linear elastic fracture mechanics (LEFM) assumes a large stress-free crack with all the nonlinearities 
of the fracture process concentrated in an infinitesimally small zone in front of the crack tip 12]. The 
basic aim of LEFM is to find the amount of energy available for crack growth and to compare it to the 
energy required to extend the crack T^, 3l|, [33] • The assumption of an infinitesimal small fracture process 



zone allows for the combination of LEFM with the framework of linear elastic mechanics of materials. 
Therefore, LEFM has been used for decades successfully for the prediction of the failure of structures. 
However, if the size of the fracture process zone is large compared to the length of the stress free crack and 
the size of the structure, LEFM results in a poor approximation of the fracture process. Most concrete 
structures are too small to fall in the range which is applicable for LEFM and a more accurate description 
by means of nonlinear fracture mechanics (NLFM) is required. NLFM describes the nonlinear response 
within the fracture process zone by means of softening stress-strain relationships. 

Constitutive models can be divided into local and nonlocal models. For local constitutive models, the 
stress at a point depends on the history of this point only. A local constitutive model for fracture of 
concrete represents an average of stochastic fracture process zones by means of a softening stress-strain 
curve. Cracks are represented by localised strain profiles within finite elements, whereby the integral of 
the energy dissipation density over the size of the finite element results in the dissipated energy. Thus, 
for local softening constitutive models for fracture, the dissipated energy depends on the size of the finite 
element. However, the energy dissipated during the fracture process of concrete is a material property and 
should be independent of the discretisation applied. This limitation of local constitutive models can be 
overcome by adjustment of the softening modulus of the softening branch of the stress-strain curve with 
respect to the finite element size [IjJ 0, [iH . With this approach, the dissipated energy is modelled mesh 
independently, as long as the inelastic strains localise in a zone of assumed size. For more complex loading 
cases, such as axial load and bending, it might happen that the inelastic strains in the compressive zone, 
which represents the fracture process zone, are independent of the element size, whereas on the tensile 
side, the width of the zones of localised strains are mesh-dependent 0- For such a case, the adjustment 
of the softening modulus does not work and more advanced methods, such as nonlocal models, need to 
be applied. Furthermore, mesh alignment can influence the orientation of strain profiles obtained [Til, lit]. 
In spite of these reservations, these models are used in many commercial finite element programs for 
modelling of fracture in concrete. 

Nonlocal models, on the other hand, are based on spatial averaging of history variables in the vicinity 
of a material point. Nonlocal constitutive models for concrete fracture describe the size of the average 
of stochastic fracture process zones independently of the size of the finite elements, as long as the finite 
elements are smaller than the width of the fracture process zone [H, 24 , 28| . This is achieved by describing 



the stress at a point of the structure by means of the weighted average of state variables in the vicinity 
of this point. Nonlocal models have several important advantages over the approach based on mesh- 
adjustment of the softening modulus. Nonlocal approaches do not rely on the assumption that strains 
localise in a mesh dependent zone. Therefore, they are capable of describing the simultaneous occurrence 
of distributed and localised failure mesh-independently. Furthermore, nonlocal models have the potential 
to describe the fracture process, in which the fracture process zone interacts with geometric boundaries, 
such as corners, stress-free cracks and topological boundaries and interfaces between different materials. 
Additionally, fracture patterns obtained with nonlocal models are not strongly influenced by the alignment 
of the finite element mesh. However, the spatial averaging of history variables is obtained by means of 
an averaging operator, which is not uniquely defined for different geometric boundary conditions [20| . 
There is generally a lack of agreement, which averaging operator should be adapted and consequently the 
nonlocal models are currently not used in commercial finite element codes, despite their advantages over 
local models relying on mesh-dependent adjustment of the softening modulus. More research is required 
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to resolve this disagreement to utilise the full potential of nonlocal models. 



Constitutive models for fracture of concrete are often based on strain-based damage mechanics, stress- 
based elasto-plasticity or combined damage-plasticity frameworks. Stress-based plasticity models are 
especially useful for the modelling of concrete subjected to triaxial stress states, since the yield surface 
corresponds at a certain stage of hardening to the strength envelope of concrete, which has a strong 
physical meaning [23, EH- Furthermore, the strain decomposition into elastic and plastic parts 



corresponds to the observed deformations in confined compression, so that unloading and path dependency 
can be accounted for in a realistic way. One of the disadvantages of the stress-based plasticity models 
is the implicit stress evaluation procedure, where for a given strain increment the final stress state must 
be iterated for satisfying the yield condition. Also, plasticity models are unable to describe the stiffness 
degradation observed in experiments. On the other hand, strain based damage mechanics models are 
based on the concept of a gradual reduction of the elastic stiffness driven by a total strain measure. 
These models are usually based on an explicit stress-evaluation algorithm, which allows for a direct 
determination of the stress state, without iterative calculation procedures. Furthermore, the elastic 
stiffness degradation observed in experiments in the case of tensile and low confined compression loading 
can be described. However, damage mechanics models are unable to describe irreversible deformations 
and are not suitable to describe the volumetric expansion ^26] . Combinations of plasticity and damage 
are usually based on isotropic plasticity combined with either isotropic or anisotropic damage [H, fl71 | . 
Anisotropic damage models for brittle materials, such as concrete, are often complex and, combined 
with plasticity, an application to structural analysis is not straightforward 0, [T3I ] . On the other hand, 
combinations of isotropic damage and plasticity are widely used and many different models have been 
proposed in the literature. One popular class of models relies on a combination of stress-based plast icity 
formulated in the effective stress space combined with a strain based damage model [2l|, [25|, |lg, Is, LLC], llJ] . 

In the present paper, both a local version based on the adjustment of softening modulus and a nonlocal 
version based on spatial averaging of strains of a combined damage-plasticity model is proposed, for a 
mesh- independent description of the fracture process of concrete. Special focus is placed on the unloading 
stiffness, which for damage-plasticity is neither perfectly secant nor perfectly elastic. To the author's 
knowledge, the mesh-dependence of the unloading stiffness for damage-plasticity models was not discussed 
before for the case of tensorial damage-plasticity models for concrete fracture. However, related work 
was presented recently for vectorial damage-plasticity models in the context of lattice approaches for 
concrete fracture [1(3] ■ The present work is not meant to favorise the use of either local models based on 
adjustment of the softening modulus or nonlocal models. For this, a much more extensive comparison 
would be required, which takes into account the influence of mesh- alignment and benchmarks in which 
multiaxial stress-states are dominant. Instead, the aim is to demonstrate how the two approaches can be 
applied to a damage-plasticity model, so that mesh-independent results are obtained. 



2 Damage-plasticity constitutive model 

The constitutive model is based on a combination of scalar damage and stress based plasticity. The 
stress-strain law is 

o- = (l-w)D e :(e-e p ) = (l-w)o- (1) 

where u> is the scalar damage parameter, ranging from (undamaged) to I (fully damaged), D c is the 
elastic stiffness tensor, e is the strain tensor, e p is the plastic strain tensor and & is the effective stress 
tensor. 

The plasticity model is based on the effective stress and thus independent of damage. The model is 
described by the yield function, the flow rule, the evolution law for the hardening variable and the 
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loading- unloading conditions. The form of the yield function is 

f(a, K )=F(a)-a y (K) (2) 
where cr y is the yield stress, which depends on the plastic hardening parameter re. The flow rule is 

e P = A- (3) 

where A is the rate of the plastic multiplier and g is the plastic potential. If the yield function coincides 
with the plastic potential (/ = g) the flow rule is associated. 

The rate of the hardening parameter re is 

dg 



|e P ]]=A 



(4) 



da 

which is the magnitude of the plastic strain rate. The yield stress is defined as 

<7y = ft (1 + H pK ) (5) 

where ft is the tensile strength. The form of H p is discussed in detail after the description of the damage 
part of the model. 

The damage part of the present damage-plasticity model is directly linked to the plasticity model. Only 
an increase of the plastic strain results in an increase of damage. However, the plasticity model is 
independent of damage. The form of the damage parameter is motivated by the definition of the inelastic 
strain for the combined damage-plasticity model. The geometrical interpretation of the inelastic strain 
and its decomposition for uniaxial tension is shown in Fig. [T] 

The stress-strain law in Eq. |T]) is rewritten as 

<t = D c : (e-(e p +u;(e-ep))) (6) 

so that the inelastic strain, which is subtracted from the total strain is 

£i = £ p +w(£-£p) (7) 

The damage history variables are chosen as the norm of the two inelastic strain components 

«di = ||e p || (8) 

and 

Kd2 — max lie — £ D || (9) 

T<t ^ 

which are the magnitude of the plastic strain and the maximum magnitude of the elastic strain, respec- 
tively. 

These history variables are related to the damage parameter for monotonic uniaxial tension, for which 
the uniaxial tension the yield equation (Eq. [2]) is written as 

f = a-a y (10) 

Thus, at the onset of yielding re<ji = and Kd2 = cr y / ' E. A linear softening law of the form 

»-*(i-S) (") 
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Figure 1: Geometrical meaning of the inelastic strain £j for the combined damage-plasticity model. The 
inelastic strain is composed of a reversible u (e — e p ) and irreversible part e p . 



was chosen, where e; is the inelastic strain in monotonic uniaxial tension. According to Eqs. CO)-©, this 
uniaxial inelastic strain is expressed by means of the two history variables n^i and n d 2 as 

Ei = K dl +UJK, d2 (12) 

Furthermore, the stress-strain law in Eq. fl} reduces in monotonic uniaxial tension to 

a = (1 - oj) Ek A2 (13) 

Setting Eq. dTTJ) equal to Eq. (O and using Eq. HU) gives 

it^di + £{En d2 - Sift 



K d2 (£{E - ft) 

This expression is further simplified by setting 



(14) 



Kd2 = a y /E = ^(l + H p K dl ) (15) 

with Eq. © for / — and n d \ = k in Eq. ([5]). This relation exists, since the history variable n d 2 is defined 
in Eq. ^ as the maximum elastic strain, which for the case of uniaxial tension is given by Eq. (|15|) . Thus, 
the damage parameter to is 

f t EK dl (1 + Ef H p ) 



(l + H pKdl )(e { Ef t -f?) 



(16) 



which depends only on K d i, which is equal to the magnitude of the plastic strain k. The hardening modulus 
Hp influences the decomposition of the inelastic strain into a permanent and a reversible component. The 
greater H p , the less is the permanent component. A physical interpretation of the hardening modulus is 
given by relating it to the ratio of permanent (plastic) and total inelastic strain, which is here denoted 
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as the unloading ratio: 

„ = (17) 

If /i = 0, the unloading stiffness is the secant stiffness. On the other hand, p = 1 corresponds to the elastic 
unloading. For uniaxial tension, the plastic strain can be determined analytically. From the consistency 
condition 

f = &-M = E(i-x)- ftHpX = (18) 

the plastic multiplier is 

A= Ei (19) 
E + f t H p 1 j 

Consequently, for monotonic loading, the total plastic strain is 

Ee ~ ft 

£ * = eTJJT v (20) 

Hence, the ratio of the plastic strain in Eq. (|2"U)) and the total inelastic strain is 

= — = Ee — — (21) 
^ ~~ £i ~ Ei (E + f t H p ) 1 ' 

For the case that the material approaches the fully damaged stage (oj 1), the total strain is equal to 
the inelastic strain, since there are no elastic strains remaining, resulting in e\ = e = Sf . Thus, Eq. (|21[) 
is written as 

»-zwtm « 22 > 

from which the hardening modulus H p is determined as 

Hp = ^d-y)-/* (23) 

For fj, = 1, which corresponds to the case of elastic unloading, the hardening modulus is H p = — 1/ef. The 
other limit, fx — > 0, corresponds to a hardening modulus H p —* oo. Furthermore, H p = corresponds to 
fi = 1—eo/sf. The present interpretation of the hardening modulus H p is based on uniaxial considerations 
and the assumption that the material is close to a fully damaged stage. For the linear softening, this 
stage is reached when s = £f . For an exponential softening law of the form 

v = f t exp(~] (24) 

which was used for the analyses of the three-point bending test, the material is not fully damaged at 
s = £f. Therefore, for exponential softening the parameter fi does not represent the ratio of permanent 
and total inelastic strain. However, the unloading ratio is independent of the element size used. For 
multiaxial stress-states the above derivations need to be extended and generalised. Some progress was 
recently made in [loj . where an approach similar to the one presented here was applied to a vectorial 
constitutive model for tension, shear and compression. 

The plastic-damage model possesses two dissipation mechanisms. Let us define the free energy as 

pip (e, £ P) J) = ^ (1 - w) (e - £ p ) : D e : (e - e p ) (25) 

Here, p is the density (specific mass), ip is the Helmholtz free energy per unit mass. Furthermore, 
isothermal processes are considered, so that the temperature remains constant and is not explicitly listed 
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Figure 2: Dissipation for (a) pure plasticity and (b) pure damage. The hatched areas denote the dissipa- 
tion for softening from A to B. 



among the state variables. The dissipation rate per unit volume is defined as 

D = a : e — pijj = tr : k 



dptp . dpip . dpip 

E+— £ n + — LU 



de " ' de p ~ p 



(26) 



The derivative of the free energy with respect to the strain results in the stress-strain relation 

-g^-e = (1 - (e - £ P ) : D : e = <t : e 
Therefore, the dissipation rate reduces to 



where the term 



and 



de p p duj 



D p = ^e p = (l~u) (£-£ p ):D c :£ p 



X> d = ^ = |(e-e t>i :D..:.r-e,,^- 



(27) 

(28) 

(29) 
(30) 



is the plastic and and damage dissipation, respectively. The two different dissipation mechanisms are 
illustrated in Fig. [2^ and b for pure damage and pure plasticity, respectively. 



2.1 Mesh-adjusted softening modulus 

The first approach to obtain mesh-independent results is based on the idea to adjust the softening 
modulus with respect to the element size. This approach, which is often called crack band approach 
0, is very common for pure damage mechanics or pure plasticity models, where the aim is to obtain a 
mesh-independent softening response. For the combined damage-plasticity model, it is also important to 
ensure a mesh-independent unloading stiffness. Therefore, both the plasticity and damage part of the 
model are made dependent of the element size. This is achieved by replacing in Eq. (|f 4p and Eq. (f2"3"]l £f 
with Wf/h, where Wf is the inelastic displacement at which the stress becomes zero and h is the size of 
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the finite element element. Consequently, the damage parameter for the crack band approach is 

/t-BKdi (l 



(l+H pKdl )(^-Ef t -f t 



(31) 



and the hardening modulus is 



h p m ( 32 ) 

For /i = la pure plastic model is obtained for which the hardening modulus is H p = —h/u>{. The smaller 
the element size h, the bigger is the hardening modulus. 



2.2 Nonlocal model 



Integral-type nonlocal models overcome the problems of the local softening models by spatial averaging 
of suitable state variables. For plasticity combined with nonlocal damage mechanics, the damage history 
variable Kdi in Eq. (|16|) is replaced by 

Kdi(x)= / " (x, a) « dl (s) ds (33) 
Jv 

so that the damage parameter u) is then determined by means of Kdi instead of Kdi- All other equations 
of the plastic damage model remain the same, i.e. no adjustment of the softening modulus with respect to 
the element size is performed. The plasticity part of the constitutive model remains local. The nonlocal 
approach is only applied to the damage part. In Eq. (|33p . a (x, s) is a nonlocal weight function that 
describes the interaction between points x and s, which decays with increasing distance between these 
points. Here, x is the point at which the nonlocal history variable is evaluated and s is a point in its 
vicinity. The weight function is usually non-negative and normalised such that 

/ a (x, s) ds = 1 for all x in V (34) 

Jv 

This can be achieved by setting 

*(x,b) = f %f I'l'iL (35) 
Jv ««, (||x-C||)dC 

Typical weighting functions are the Gauss-like function or the bell-shaped truncated polynomial 
function. The bell-shaped function 

« h J R (r)=\ K 1 "^) 2 [i ^ R (36) 
if r > R 

which has a bounded support, i.e. it vanishes at distances r = ||x — s|| exceeding the interaction radius R. 
The scaling factor is c = 16i?/15 and c = irR 2 /3 for the one and two-dimensional case, respectively [32j |. 
The present combination of local plasticity with nonlocal damage is computationally more efficient than 
standard integral-type nonlocal plasticity models, since the plasticity part, which is implicit, remains 
local and only the equivalent strain in the explicit damage part is averag ed [H, M, Q. In [H, Q it 
was concluded that mesh-independent results can only be obtained, if the hardening modulus H p of the 
plasticity model is greater than the critical hardening modulus of the material. Therefore, only /i-values 
corresponding to H p > are considered in the present study. 
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Figure 3: Constitutive response: (a) Stress-strain curves for /i = 0.01, /i = 0.5 and /i = 0.99, (b) Dissi- 
pation versus strain for fi = 0.01, /i = 0.5 and [i = 0.99. 
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Figure 4: Geometry of the uniaxial tensile specimen. 



3 Constitutive response 



For illustration of the above concepts, the interaction of damage and plasticity is here studied for the 
case of softening in uniaxial tension. The material properties are chosen as E = 20 GPa, / t = 2 MPa 
and £f = 0.001. In Fig. [3^,, the stress-strain relations are shown for fi — 0.99, 0.5 and 0.01. For [i = 0.99 

(Hp « ), the damage variable u> in Eq. (|16p is equal to zero and a plastic response (e; w e p ) is 

obtained. On the other hand, for \x — > (H p — > oo) the plastic multiplier in Eq. (|19p is zero and a 
damage response with £; « ujs is obtained. The evolution of the dissipation determined from Eq. (|28[) 
for the three values of /i is shown in Figure^. For the linear softening law in Eq. (J5J), the dissipation 
evolution of the plasticity model is nonlinear, whereas for the damage model a linear dissipation evolution 
is obtained. The dissipated energy per unit volume for the present linear softening law is determined to 

9f = 7r/t£f- 



4 One-dimensional damage-plasticity response 

Two ways of resolving the issue of mesh dependency associated with local softening constitutive models, 
namely mesh-adjusted softening modulus and nonlocal damage models, are discussed further. The two 
approaches are applied to a bar subjected to uniaxial tension shown in Figured The material parameters 
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Figure 5: Effect of mesh refinement on the numerical results: a) load-displacement diagrams, b) strain 
profiles. 




Figure 6: Nonlocal plastic-damage model: (a) Load-displacement curves for 21, 41 and 81 elements and 
b) evolution of the strain distribution for 81 elements at four stages of analysis. 



are the same as used in the previous sections. The load-displacement curves for three mesh-sizes with 1, 
5 and 21 elements are shown in Figure EK- The element located in the middle of the bar of the meshes 
is weakened by 0.1%. The strain profiles are shown in Figure Od. The load-displacement curves are 
independent of the mesh-size. However, the strain profiles depend on the mesh size. 

The nonlocal model is applied to the same uniaxial problem which was analysed in the previous section 
with the mesh-adjusted softening modulus. Three different meshes with 21, 41 and 81 elements were used. 
A section of 0.005 m in the middle of the bar was weakened (Figure [H]) by 0.1%. The nonlocal radius 
was chosen as R — 0.01, which is approximately twice the element size of the coarsest mesh. The linear 
softening curve from Equation (0) with e f = 0.002 was used. The other material parameters were the 
same as in the section presenting the constitutive response. The unloading ratio was chosen to /i = 0.5. 

The load-displacement curves are shown in Figure [5^ and demonstrate an almost mesh-independent 
result. Furthermore, the strain profiles in Figure [7] are almost identical for the three meshes. However, 
the load-displacement curves are nonlinear, although a linear local softening curve was used in Eq. (llip . 
The nonlinear load-displacement response is explained by the shrinking of the zone of localised strains 
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Figure 7: Nonlocal plastic-damage model: a) Strain distribution for 21,41 and 81 elements at stage 3 
marked in Fig. [5^. 




(Fig. [5b)- With decreasing load, the strain localises further, so that the load-displacement curve becomes 
steeper. Furthermore, the influence of the ratio on the fracture process zone and load-displacement curves 
is investigated. For the nonlocal model, the damage part is nonlocal and the plasticity is local. To ensure 
mesh-independent results, the hardening modulus H p must be greater than zero. This corresponds for the 
chosen model parameters to the condition fi > 0.9. Four unloading ratios of = 0.01, 0.5, 0.75, 0.9 were 
used for the analysis with the fine mesh. The load-displacement curves and strain profiles are shown in 
Fig. [51 The softening branch of the load-displacement curves differs slightly. However, the strain profiles 
in Fig. [SJd are almost independent of the parameter \i. 



5 Three point-bending test 

The softening modulus adjustment and nonlocal approach are here applied to the analysis of the three 
point bending test, for which the experimental results were reported in (23[. For this analysis, the two- 
dimensional Rankine yield surface with an associated flow rule was used for the plasticity part of the 
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Figure 9: Geometry and setup of the three-point bending test. 



damage-plasticity model |32| ■ The hardening law of the plasticity part was chosen according to Eq. ([2"3"jl . 
For the damage part, an exponential softening law was chosen. The geometry and the setup of the three 
point bending test is shown in Fig. [5] Three finite clement meshes with element size 7.5, 5 and 2.5 mm 
in the central part of the specimen were used. The coarse mesh is shown in Fig. [5J The three point 
bending test was analysed with the damage-plasticity model using the mesh-adjusted softening modulus 
and the nonlocal model. The model parameters for the mesh-adjusted softening modulus were chosen 
as E = 20 GPa, / t = 2.4 MPa, Wf = 0.000047 m. Furthermore, the Poisson's ratio required for the 
two-dimensional analysis was chosen as v = 0.2. For the nonlocal approach the critical crack opening Wf 
was replaced by the critical strain Ef — 0.001025 and the nonlocal radius is set to R = 0.025 m. These 
parameters resulted in satisfactory results for a damage-plasticity approach reported in [7f, in which 
damage was determined by the plastic strain and the plasticity model was chosen to be perfect plastic 
(H p = 0), which corresponds to \i — 1 — ft/ (Eef) — 0.883. In Fig.fTUlthe present two modelling approaches 
with fi = 0.883 are compared to the experimental results reported in [23| . The results obtained with the 
two modelling approaches are in satisfactory agreement with the experimental results. In the following 
sections the mesh- dependence and the influence of the unloading ratio /i arc investigated. For the mesh- 
adjusted softening modulus, the load-displacement curves for the three meshes are shown in Fig.fTTa. The 
results are almost independent of the size of the elements. Furthermore, the load-displacement curves 
for the ratios /j, — 0.01, 0.5 and 0.9 are presented for the fine mesh in Fig. [12b . The unloading part of 
the load-displacement curves differs depending on the chosen unloading ratio //. However, the softening 
part is almost independent of this ratio. For the nonlocal model, the load-displacement curves for the 
three meshes are shown in Fig. [I]J]for /i = 0.5. Again, the load-displacement curves are independent of 
the element size. However, the maximum load is considerably higher than for fj, = 0.883, for which the 
nonlocal model was calibrated to match the experimental results. This overestimation of the peak load is 
further investigated by varying the unloading ratio (i. The load-displacement curves for /i = 0.5, 0.8 and 
0.883 are shown in Fig. [TU The smaller the fi, the greater is the peak load. This strong dependence 
on the unloading ratio /i is caused by the presence of the notch in the problem considered. The nonlocal 
history variable of the damage model, which is used to determine the damage parameter, is computed by 
means of the spatial average of the local history variables in the vicinity of the point for which the stress 
is evaluated. The peak load of the beam is very sensitive to the stress-strain response of the material 
directly above the notch. For the stress-evaluation of this region, the areas below the notch are included 
in the averaging. However, these areas are not strained during the loading process, so that the local 
history variable of the damage model will remain zero in these regions. Therefore, the material above the 
notch is stronger than points further away from the notch along the ligament. For fi — > 0.883 the stress 



12 



2.5 






0.4 0.6 
d [mm] 



Figure 11: Three point bending test: Load-displacement curves obtained with the mesh-adjusted softening 
modulus for three mesh sizes with \i = 0.5. 
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Figure 12: Three point bending test: Load-displacement curves obtained with the mesh-adjusted softening 
modulus for three ratios of /j, = 0.01, 0.5 and 0.99 for the fine mesh. 
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Figure 13: Three point bending test: Load-displacement curves obtained with the nonlocal model for 
three mesh sizes with /j, = 0.5. 
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Figure 14: Three point bending test: Load-displacement curves obtained with the nonlocal model for 
three ratios of [i = 0.5, 0.8 and 0.883 for the fine mesh. 



is limited by the local plasticity model. On the other hand, for lower values of /i, the plasticity model 
exhibits hardening. Thus, a reduction of the unloading ratio results in an increase of the strength of the 
region above the notch and, consequently, in an increase of the peak load of the beam. 



6 Conclusion 

The present work deals with a damage-plasticity approach for the analysis of fracture in concrete. Two 
approaches to obtain mesh-independent load-displacement curves were applied to a uniaxial bar subjected 
to tension and a three-point-bending test. The first approach relies on an adjustment of the softening 
modulus with respect to the finite element size, whereas the second approach is based on nonlocal aver- 
aging of damage history variables. The following conclusions can be drawn: 

• For both approaches the present damage-plasticity model results in a mesh independent description 
of the softening and unloading regime of the load-displacement curves. 

• For the adjustment of the softening modulus, a change of the unloading stiffness does not influence 
the softening regime of the load-displacement curves for the finite element meshes used in the 
present study. 

• For the nonlocal model, the choice of the unloading stiffness influences strongly the load-displacement 
curves obtained for the three-point bending test. This dependency is explained by the influence of 
the boundary near the notch, which affects the nonlocal averaging procedure. 

The present results indicate that the nonlocal damage-plasticity model exhibits a strong boundary effect. 
More research is required to develop boundary operators that overcome this influence of boundaries on 
nonlocal averaging. This is one of the reasons why nonlocal modes are currently not yet implemented 
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in commercial finite element codes. Nevertheless, nonlocal models are more general than models based 
on the adjustment of the softening modulus with respect to the element size. They are less sensitive to 
the alignment of the finite element mesh and produce even mesh-independent results, if localised and 
distributed failure modes coexist Ml. 
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